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WATER  ENTRY 

This  report  describes  a method  for  predicting  pressures,  forces,  and  moments 
on  arbitrary  bodies  during  vertical  and  oblique  water  entry.  Also  Included 
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INTRODUCTION 

A common  problem  in  the  design  of  bodies  which  enter  the  water  at  high  speeds 
is  the  determination  of  the  surface  pressures,  forces  and  moments  during  water 
impact.  This  paper  describes  an  engineering  method  for  calculating  these 
quantities.  A simplified  potential  model  is  used  which  replaces  the  water's 
free  surface  with  an  effective  planar  surface  that  is  positioned  using  un 
empirical  parameter  available  in  the  literature  for  a wide  variety  of  shapes. 

To  confirm  predictions,  calculations  are  compared  to  experiment  for  the 
oblique  entry  of  spheres,  disks,  and  ogives  and  for  the  vertical  entry  of 
spheres,  cones,  ogives,  and  cusps.  Surface  pressures  agree  well  with  measurement 
reflecting  both  the  model  geometry  and  location  on  the  model.  The  calculated  drag 
and  lift  exhibit  close  agreement  with  experimental  vulues,  particularly  prior  to  the 
peak  loads.  At  later  times  the  shape  of  the  hydraulic  cavity  muBt  be  taken 
into  account  and  an  approximate  procedure  for  doing  this  is  described.  A 
program  listing  and  sample  runs  are  provided  hs  well  as  instructions  for 
using  the  code. 

Attempts  to  analyze  the  water-entry  problem  originate  circa  1929  with  the 
work  of  von  Karman1.  Comprehensive  surveys  of  this  field  are  provided  by  May2, 
Thigpen5,  Szebehely^,  and  Moran5,  The  main  thrust  of  early  work  follows  the 
formulation  developed  by  von  Karman  and  Wagner^,  In  this  approach  a potential 


■*von  Karman,  T. , "The  Impact  on  Seaplane  Floats  During  Landing,"  NACA  TN  321, 

Oct  1929 

2May,  Albert,  "Forces  at  Water  Impact,"  Alden  Research  Laboratories,  ARL  119-72/SP, 
Dec  1972 

^Thigpen , A.,  "Water-Entry  Technology  - A Review,"  Sandia  Corporation  Technical 
Report  SC-Dr  71  0196  (Jun  1971) 

^Szebehely,  V.  G.,  "Hydrodynamic  Impact,"  Appl„  Mech,  Rev.,  12,  297-300,  1959 

5Moran,  J.  P,,  "On  the  Hydrodynamic  Theory  of  Water  Exit  and  Entry,"  Therm 
Advanced  Research  Technical  Report  TAR-TR  6501  (Mar  1965) 

^Wagner,  H.,  "liber  Stoss-und  Gleitvorgange  an  dcr  Oberflache  von  Flussigkelten," 
ZAMM  1_2,  4,  193-215,  1932 
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flow  modal  Is  used  and  forces  are  calculated  by  the  added  mass  concept.  The 
submerged  portion  of  the  body  Is  often  ftcued  or  replaced  by  another  with  the 
same  surface  cross-sectional  area  tor  which  a closed  form  solution  is  available. 

A linearlized  version  of  the  free-surface  boundary  conditions  is  applied  to 
determine  the  surface  shape.  Most  of  the  theories  are  restricted  to  vertical 
entry  of  simple  geometries.  In  recent  years,  computational  efforts  have  been 
made  to  obtain  a solution  using  the  non-linear  boundary  conditions.  An  early 
example  of  such  work  is  that  of  Chu  and  Falconer7.  A relaxation  method  was  used 
to  solve  the  potential  problem  for  arbitrary  bodies.  This  project  was  abandoned 
due  to  problems  with  excessive  computational  time  and  surface  contact  discontinuities 
The  same  formulation  for  the  vertical  entry  of  cones  has  been  treated  by  Weber® 
using  a distribution  of  source  dipoles.  More  recently,  Shere  and  Vender  Vorst9 
and  Vender  Vorst  and  Rogers^  have  used  the  marker  and  cell  method  to  develop  a 
detailed  viscous  model  of  vertical  cone  entry. 

The  objective  of  the  current  study  has  been  to  develop  an  engineering  tool 
for  calculating  pressures  and  loads  which  is  simple  to  use,  inexpensive  to 
exercise,  and  applicable  to  a wide  variety  of  geometries.  Accordingly,  the 
philosophy  of  the  current  program  has  been  to  include  only  those  portions  of  the 
problem  which  can  be  shown  empirically  to  be  necessary.  The  current  work  combines 
a simple  flow  field  model  with  the  potential  flow  computational  techniques  of  Hess 
and  Smiths  to  form  an  extremely  versatile  approach  which  can  be  applied  to  a 
wide  variety  of  geometries  over  a broad  range  of  entry  conditions.  The  success 
of  this  calculative  method  indicates  that  a detailed  description  of  the  free 
surface  is  not  necessary  for  the  purposes  of  calculating  entry  pressures  and 
loads. 
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PROBLEM  FORMULATION 

The  flow  field  about  the  entering  body  is  assumed  to  be  described  by  a 
potential  model.  The  free  surface  is  simulated  with  an  effective  planar  surface 
whose  location  is  defined  using  the  measured  wetting  factor  C„.  This  parameter  is 
equal  to  the  ratio  of  h/h'  where  h is  the  effective  depth  of  the  model  and  h*  is 
its  actual  depth  or* penetration  below  the  original  free  surface  (see  Fig„  1).  The 
governing  equations  and  boundary  conditions  are: 


Governing  equation: 
Boundary  conditions: 

7^4  “ 0 

(1) 

(a)  On  body  surface: 

— V<b« e - V „ • e 
n E n 

(2) 

(b)  On  the  effective  planar  Burface: 

1 

( V - -(C  -1)V  » k 

| 8 W p 

| 

(3a) 

1 

' 4-0 

(3b) 

Pressures  are  calculated  from  successive  solutions  at  differing  depthB  using 
the  unsteady  Bernoulli  equation: 

s-jp*  - - \ m)2  (m 

This  equation  must  be  case  in  body  fixed  coordinates  since  4 is  calculated  at  the 
same  point  on  the  body  in  successive  steps.  Thus,  equation  (4)  becomes: 

V*  ■ If  - v*  - <5) 

The  above  problem  can  be  put  in  nondimens ional  form  by  applying  the  following 
transformations: 

A 

4 - 4/VjD 

AAA 

x,  y,  z - x/D,  y/D,  z/D  (6) 


A A 

h - (V  »k)C  tV./D 
p w I 
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Now  equations  (1)  through  (3)  become: 


V <p  = 0 


-(V  $)*C 


e " V_ 
n E 


V - -(C  -1)V  “k 
s w p 


* - 0 


(7) 

(8) 


on  effective  planar  (9) 
surface 


The  nondimens ional  pressure  is: 


C - 2 -gr  [<V  *k)tC  32*  - 2V  "7l  - 0?J)2 

p 3t  p w 9h  E 

For  constant  entry  conditions  (i.e,,  V and  C are  fixed)  the  above  becomes: 

p w 


C “ 2C  sin8  - 2 e ‘Vcp  - (Viji)*1 


(10a) 


(10b) 


Dh 


These  two  equations  indicate  that  the  calculated  pressure  and  force  coefficients 
are  independent  of  the  model  and  entry  velocity  scale  (i„e„,  D and  Vj  respectively) „ 
The  value  of  these  two  parameters  must  be  simulated  through  an  appropriate  choice 
of  C „ Also  it  is  evident  that  for  constant  entry  conditions,  depth,  not  time  is 
the  most  natural  independent  variable. 


The  boundary  conditions  used  in  the  current  study  are  similar  to  the  linearized 
version  applicable  to  slender  bodies.  The  linearized  conditions  are  that: 


<j>  ■ 0 on  z ■ 0 


(11) 


V 

B 


- H (x,y,0)  k 


(12) 


These  conditions  follow  from  the  nonlinearized  form  by  dropping  the  quadratic 
terms  which  are  second  order  as  long  as  <p  and  its  derivatives  are  small  near 
the  surface.  The  present  model  applies  an  empirical  correction  to  the  surface 
velocity  described  by  equation  (12), 

POTENTIAL  FLOW  SOLUTION 


At  each  depth  the  problem  requiring  solution  is  described  by  equations  (7) 
through  (9)  and  is  directly  amendable  to  the  potential  flow  techniques  developed 
by  Hess  and  Smith  which  use  a distribution  of  sources  and  sinks.  The  surface  of 
the  body  under  consideration  is  divided  into  quadrilaterals  and  a constant 
source  strength  is  assumed  to  exist  throughout  each  element.  The  source 
strengths  are  determined  by  satisfying  equation  (8)  at  the  centroid  of  each 
element  which  results  in  a system  of  N simultaneous  equations  of  the  form: 


03) 


11 
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Here  oj  Is  the  source  strength  of  element  j,  e is  the  unit  vector  normal  to 

n 1 

element  l,  and  A.,  Is- the  normal  velocity  induced  on  element  i by  unit  source 
strength  on  element  j.  Equation  (13)  is  solved  directly  using  the  method  of 
reference  12,  When  the  number  of  elements  exceeds  120,  solution  is  accomplished 
using  a series  of  blocks. 


Once  the  source  strengths  are  determined,  the  velocity  and  potential  at  the 
centroid  of  each  element  can  be  calculated: 


V 

n 


i 


<t>  - 


oj 


(14) 


V 


Yl 


■ e 

v 


e 

n 


i 


Here  B^.p  C . . , and  D,,  represent  the  quantities  V.,  V and  <J>  induced  on  element  i 

by  element  Jjassumingj element  j has  a source  strerigthnof  one.  The  term  of 

matrices  [A],  [B],  [C],  and  [ D]  are  evaluated  using  the  closed-form  expressions 

given  in  reference  11  which  are  reproduced  in  Appendix  A.  Equations  (13)  and  (14) 

are  cast  in  the  inertial  frame  of  reference  where  V » 0, 

00 

In  applying  the  above  method  to  the  water-entry  problem,  only  the  submerged 
portion  of  the  body  (i.e.,  below  the  effective  planar  surface)  is  considered. 

The  extra  condition,  tj>  * 0,  is  satisfied  on  the  effective  planar  surface  by 
locating  image  elements  above  this  surface  as  shown  in  Figure  2„  The  strength  of 
the  image  element  is  equal  in  magnitude  but  opposite  in  sign  to  the  original  one,. 


12. 


Forsythe,  G.,  and  Moler,  C.,  Computer  Solution  of  Linear  Algebraic  Systems, 
Prentice-Hall,  Englewood  Cliffs,  NJ,  1967 
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If  flic  entry  body  possess  symmetry  about  the  v'-z'  plane,  on  1 v half  of  the  model 
is  gridded  since  symmetric  element  pairs  have  the  same  source  strength.  For 
such  a body,  four  types  of  elements  have  the  same  source  strength  magnitudes  and 
their  influence  coefficients  are  grouped  together.  The  terms  A^. , . t C^,  and 

D reflect  the  influence  on  element  i of  element  j , its  linage,  the  corresponding 
symmetric  element  and  its  image.  If  the  entry  body  does  not  possess  planar 
symmetry  the  entire  face  must  be  gridded.  Here  each  influence  coefficient 
reflects  only  the  effect  of  an  element  and  its  image. 

COMPUTATIONAL  PROCEDURE 

A series  of  points  or  nodes  are  defined  on  the  surface  of  the  body  of 
interest  in  x*,  y',  z*  coordinates.  These  are  arranged  into  groups  of  four  to 
form  planar  quadrilateral  elements  as  shown  in  Figure  3,  The  several  different 
options  available  for  defining  nodes  and  elements  on  arbitrary  bodies  are 
discussed  in  Appendix  C. 

Solution  of  the  Potential  Problem 


The  entire  body  may  have  an  nribtrary  entry  velocity  and  rotation  in  the  y-z 
pLane  is  allowed.  The  computation  proceeds  by  inserting  the  model  into  the  water  in 
a series  of  steps  each  at  n depth  greater  than  the  previous  one.  The  entry 
velocity  and  the  increment  in  model  depth  can  be  varied  from  step  to  step.  At 
every  step  the  group  of  elements  comprising  the  submerged  portion  of  the  model 
are  redefined  and  arranged  into  a form  amenable  to  the  ealeulative  procedure 
outlined  in  the  previous  section.  The  nodes  definining  a particular  element  are 
checked  to  determine  whether  they  are  above  or  below  the  water  line.  Elements 
with  nil  four  nodes  above  the  water  line  are  discarded  while  those  with  all  four 
below  it  are  Included  without  change.  Element  which  are  intersected  by  the  water 
surface  may  have  either  one,  two  or  three  submerged  nodes  as  shown  in  Figure  4,. 

In  all  cases  two  new  nodes  are  generated.  Given  two  nodes,  one  below  the  water 
surface  (x^,  y ^ , Zj)  and  one  above  it  (X2,  y^,  t'":>  nt’w  nocl1'  located  at  the 

water  surface  on  <1  line  intersecting  these  two  points  is: 


z - 0 
new 


ynew  " yl  + (y2“yl)  T*7^) 


X “ X.  + (x--x,)  r~ \" 

new  ! '21  ( z - z2> 


(15) 


When  an  element  has  only  one  node  submerged,  it  is  necessary  to  define  a third 
new  node  in  order  to  obtain  a quadrilateral  element.  This  last  new  node  Is  placed 
midway  along  the  surface  edge  of  the  element.  If  only  one  node  Is  above  the 
water  surfc.ce  the  generation  of  two  new  nodes  resu.l  ts  in  a pent  a 1 nt era  1 element,. 

Here  again,  a third  new  node  is  added  and  the  element  Is  broken  into  two  parts  each 
of  which  is  now  quadrilateral. 
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It  is  necessary  to  define  a set  of  element  coordinates  associated  with  each 
quadrilateral  element  used  in  the  computations.  This  h,  £,  Y coordinate  system 
is  shown  in  Figure  5 and  the  corresponding  unit  vectors  are  as  follows: 


i 


i 


I 

(i 


r 

i 

i 


I 


(x3-x1)i  + <y3-y1)j  + (z3-z1)k 
^ y<x3“x1)2  + ^y^y].^  + (z3-zi^ 

i-  X[ (x  -X4)I  + (y2-y4)j  + (z2-z4)k]  (16) 

e * -*» — - ■ ■ ■ ■ ■ — _ ~ — 

Y |e^x[ <x2-x4)i  +- (y2-y4)j  + (z2-z4)k]| 

e"  m e"  Xe\. 
n r K 

Here  the  subscripts  refer  to  the  corner  numbers  shown  in  Figure  5, 

At  every  step,  equations  (1)  to  (3)  are  aolved  using  the  potential  flow 
method  discussed  in  the  last  section.  The  value  of  the  velocity  and  potential 
at  each  element  centroid  is  stored  for  future  use  in  determining  C . In  the  case  of 
elements  which  have  been  split  into  two  (Fig.  4c),  a single  area  weighted  average 
value  is  retained. 

Calculation  of  Surface  Pressures 


At  each  depth  the  pressure  coefficient,  C , is  evaluated  at  each  element 
centroid  using  equation  (5)  which  is  in  a bodj?  fixed  frame  of  reference: 


C 

P 


P-Po= 


2 36  . 2VV 

V 2 V 2 

I VI 


(17) 


The  fluid  velocity, gV,  which  appears  in  this  equation  is  directly  determined 
at  each  depth,  but  muBt  be  calculated  using  the  value  of  cp  at  the  same 
body  locations  in  adjacent  steps.  The  general  expression  used  to  calculate 
this  quantity  at  the  n^h  step  is: 


if 


n 


At 


n-1 


Vi + <Vi  - ♦«-!>  Tat^-; 


n+l  + At  . ] 
n-i 


ij>r  - <pr  - 1 

cn  cn 


•n-l 


At 


n-1 


(18) 
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n-*-l 


Kcn+1 


At 


n+l 
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Here  i|>  is  the  value  of  the  potential  at  the  element  centroid  where  the  pressure 

n (.ji 

ia  being  calculated  at  the  n step.  The  quantity  At  , is  the  time  interval 
between  steps  n-1  and  n.  Similarly,  At  is  the  time  interval  between  steps  n 
and  n+1.  Note  that  if  At  , ■ At  . 1 , tReabove  expression  reduces  to  the 
central  difference,  n n 


4» 


n+l 


cn-l 


2At 


n+l 


(19) 


94> 


n for  elements  which  are  modified 


Special  problems  arise  in  calculating  7t 
(i,,e,,  intersected  by  the  water  surface)  in  any  of  steps  n-1,  n,  and  n+l.  This 
is  because  the  body  fixed  coordinate  of  a modified  element  centroid  differs 
from  its  unmodified  value  and  hence  4>  Is  not  known  at  the  same  point  on  the 
body  surface  for  the  required  number  of  adjacent  steps.  To  handle  this  situation 
local  similarity  is  assumed.  ThiB  assumption  holds  that  at  any  point  within  un 
element  $ is  only  a function  of  the  length  of  time  that  this  point  has  been 
submerged..  This  removes  the  necessity  of  knowing  <j>  at  the  same  point  on  the 
body  surface*  Hence  the  values  of  associated  with  the  same  element  centroid 
are  used  in  equation  (18)  regardless  of  whether  the  element  is  modified  in  any 
of  the  three  required  adjacent  steps,.  If  an  element  is  modified,  the  associated 
time  interval  between  it  and  proceeding  or  following  steps  to  be  used  in 
equation  (18)  is  no  longer  the  time  interval  between  successive  steps.,  The 
required  time  interval  to  be  used  in  place  of  Atn_^  is: 


At  - 


h - h . 
n n-1 

C'  V 
w z 


(20) 


An  analogous  expression  applies  for  determining  At  +.,  Here  h is  the  depth  of 
the  element  centroid  at  step  n while  C'w  and  arc  the  wetting  factor 

and  the  z velocity  component  of  the  element  centroid  between  steps  n-1  and  n„ 
Since  the  model  may  rotate  In  the  y-z  plane  the  velocity  vector  of  different 
points  on  the  body  surface  will  vary  with  location,.  Hence,  the  wetting  fnrtor 
C'  used  In  equation  (20)  is  a local  value  and  not  that  prescribed  Tor  the 
en¥ry  body,.  Between  steps  n-1  and  n this  parameter  is  calculated  from: 


C 


i 

w 


V + (V'k) (C-l) 
z p w 

V 


(21) 


where  C is  the  prescribed  time  interval  and  wetting  factors  between  n-1  and 
w 

n.  An  analogous  expression  is  used  to  determine  C between  steps  n and  n + 1. 


lr) 
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For  modified  elements  located  near  the  water  surface  it  is  also  possible 
to  use  the  boundary  condition  <}>  . ■ 0 at  h ■ 0.  This  condition  must  be 

applied  if  the  element  at  whichntne  pressureiB  to  be  calculated  is  not 
present  in  the  preceding  step  (i„e.,  no  part  of  it  was  submerged).  It  has 
also  been  found  advantageous  to  use  this  condition  for  elements  modified  in 
step  n. 

The  local  similarity  assumption  is  strictly  applicable  for  the  oblique 
entry  under  constant  velocity  and  orientation  of  an  infinite  plate.  For  plates 
of  infinite  length  but  finite  and  constant  cross-sectional  geometry  this 
assumption  holds  in  x ■ constant  planes  shown  in  Figure  6.  This  assumption 
is  well  founded  for  bodias  where  conical  similarity  is  applicable  if  for 
successive  steps  Ah  « h.  On  three-dimensional  models  this  assumption  is  most 
accurate  on  portions  of  the  model  where  the  surface  geometry  varies  slowly. 


The  pressure  coefficient  on  the  model  at  the  water  surface  is  singular. 
This  can  be  seen  by  casting  equation  (4)  in  a frame  of  reference  moving  with 
the  effective  planar  surface. 


C 


■ 

_2_  M+  v 
2 It  a v 2 

VT  L J 


3d) 

On  this  surface  ■ 0.  Due  to  a source  discontinuity  at  the  intersection  of 
the  model  and  the  water  surface  V - "“and  C Fortunately,  the  value 

of  C recovers  quickly  with  depth  and  assumed  a positive  value  well  before 
the  experimentally  observed  pressure  peak.  For  the  element  sizes  uBed  in 
this  study  the  first  value  of  pressure  calculated  for  each  element  1b 
usually  positive. 


Negative  C values  can  also  be  obtained  on  the  sides  of  the  entering  body 
if  allowance  ii  not  made  for  the  flow  cavity.  Such  values  are  set  to  zero 
for  the  purposes  of  calculating  total  model  loads. 


USE  OF  THE  NUMERICAL  MODEL 


In  using  the  described  numerical  model  it  is  necessary  to  specify  the 

wetting  factor,  Cw,  the  increment  in  depth  between  successive  steps  (Ah)  and 

to  construct  an  appropriate  grid.  In  some  cases  it  is  also  advisable  to  apply 

a correction  to  the  pressure  calculated  on  modified  elements.  The  parameter 

Cw  can  be  determined  from  a body  of  existing  experimental  data  which  will 

be  reviewed  in  conjunction  with  specific  applications.  The  numerical  effect 

on  C_  and  C of  varying  C is  illustrated  in  Figures  7 and  8 respectively, 
k'oo  P W 

With  decreasing  values  of  is  reduced  in  magnitude  and  peak  values 

00 

occur  at  a later  time.  The  peak  pressure  coefficient  Increases  with  increasing 

values  of  C , but  its  rate  of  decay  is  also  accelerated, 
w 
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Selection  of  Ah 


The  flow  field  properties  at  any  particular  depth  are  Independent  of  solutions 
at  other  depths  and  hence  of  Ah.  However,  calculation  of  C , as  discussed  In 
the  previous  section,  requires  values  of  4>  from  adjacent  stlps.  In  as  much  as 
the  present  method  calculates  only  a single  pressure  for  each  element,  it  is 
desirable  that  this  pressure  represent  an  average  for  the  element.  To  ascertain 
an  appropriate  step  size  for  accomplishing  this,  the  constant  velocity  and 
orientation  entry  of  the  flat  plate  of  finite  length  shown  in  Figure  6 is 
considered.  Defining  F?  to  be  the  depth  of  an  arbitrary  point  p on  the  surface 
of  the  plates 


c»*ln9v-  I [ 


Transforming  the  above  into  the  x",  y",  z"  coordinate  system  of  Figure  6 which 
is  fixed  on  the  effective  planar  surface: 


H ’ " (CwSil0V  [■3yrr  p 


The  average  value  of  -g—  for  the  rectangular  elements  shown  in  Figure  6 with  a 
paii  of  edges  parallel  to  the  water  surface  Ibs 

V " v " 

2 y2 

C sinOV  4 ' 


1 - -&1  fj  „ If  *■'*"  - - *,(l,i>1 

xl"  yl"  2 1 


(23) 


or 


30  " »*(tl> 


Tt  ‘ 


(t2  - tx) 


Here  and  cp* ^ 2 ) are  t*ie  avera8e  values  of  0 along  the  upper  and  lower 

edges  of  the  element  respectively , A is  the  element  area,  and  t„  and  t.  are  the 
lengths  of  time  that  the  lower  and  upper  edges  of  the  element  have  been  submerged. 

Equation  (23)  is  of  the  same  form  as  the  central  difference  expression 

used  in  evaluating  (equ.  (19))u  If  the  pressure  on  the  rectangular  element 
of  Figure  6 is  being  evaluated  at  step  n,  equations  (19)  and  (23)  become 
Identical  if  the  step  size  is  chosen  such  that  at  step  n-1  and  n+1  the  element 
centroid  Lies  at  the  top  and  bottom  edges  of  the  element  In  step  n respectively,, 
Hence  the  step  size  should  be  chosen  ho  that  each  element  is  completely  submerged 
In  two  steps  or 

Ah  - <y 2 1 - y11)cos0/2 
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The  preceedlng  analysis  is  not  strictly  applicable  to  three-dimensional 
bodies  entering  with  variable  velocity  and  orientation  and  composed  quadri- 
lateral elements  with  edges  not  necessarily  parallel  to  the  water  surface. 

However,  it  is  taken  us  a guide  and  Ah  is  picked  to  insure  that  the  average 
element  is  submerged  in  two  steps.  Thus, 

Ah  ~ S-  cos  G * /2  (24) 

where  & is  the  element  characteristic  length  and  0'  is  the  typical  element 

orientation  angle  with  respect  to  the  vertical.  This  criteria  is  easily 

applied  on  flat  plates  or  cones.  Spheres  and  other  bodies  with  curvature  in 

the  axial  direction  are  more  difficult  to  deal  with.  The  above  1b  satisfied 

only  approximately  with  elements  perpendicular  to  the  direction  of  motion 

being  most  heavily  weighted.  In  practice,  it  is  these  elements  which 

experience  the  largest  pressures  and  hence  are  most  crucial  to  the  problem  solution. 

On  bodies  with  curvature  in  the  axial  direction  the  size  of  the  elements  1b 

increased  as  their  orientation  approaches  the  direction  of  motion. 

The  effect  of  varying  grid  size  and  hen^c  Ah  on  Cj^  and  C is  illustrated 
in  Figures  9 and  10  respectively.  On  flat  vurl.ices  an  accurat!  solution  is 
obtained  with  only  a email  number  of  elements  Thr  principal  effect  of 
increasing  the  number  of  elements  is  to  reduce  the  peak  calculated  pressure. 

Since  these  peak  pressures  act  over  small  areas  the  drag  coefficient  Is 
relatively  insensitive  to  increases  in  the  number  of  elements.  More  complex 
shapes  naturally  require  the  use  of  a larger  number  of  elements. 

Describing  the  Entry  Body  with  Quadrilateral  Elements 

In  setting  up  a grid,  beat  results  are  obtained  if  the  afterbody  is  neglected 
and  only  the  nose  of  the  entering  body  is  gridded.  The  pressures  and  source 
strengths  on  afterbody  elements  are  small.  On  models  with  well  rounded  shoulders 
exclusion  of  these  elements  decreases  the  required  computational  effort  without 
strongly  effecting  the  solution.  On  models  with  sharp  edges  such  as  disk 
cylinders,  inclusion  of  the  afterbody  elements  imposes  the  requirement  that  the 
flow  make  a sharp  turn  about  the  edge  of  the  face.  This  requirement  1h  physically 
unrealistic  since  the  flow  will  separate  at  thu  face  edge.  Use  of  afterbody 
elements  in  this  case  increases  the  flow  velocity  on  elements  near  the  edge 
which  in  turn  decreases  the  calculated  pressures..  This  is  illustrated  in 
Figures  11  and  12  which  give  the  calculated  drag  for  the  oblique  entry  of  a 
disk  cylinder  w.th  ind  without  afterbody  element  respectively.  Results 
obtained  in  Figure  12  without  the  use  of  afterbody  elements  are  in  much  closer 
agreement  with  experimental  results  and  required  a smaller  amount  of  computational 
effort. 

Cavity  Modeling  Using  No  Load  Elements 

At  time  following  peak  impact  loading  the  existence  of  the  flow  separation 
region  or  cavity  about  the  afterbody  of  entry  models  may  have  a significant  effect 
on  the  model  surface  pressures  and  positive  steps  must  be  taken  to  model  it. 
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Accordingly,  no  load  elements  have  boon  lnt roduood  Into  the  computations. 

Those  elements  are  placed  on  the  water-cavity  Interface  and  their  purpose 
Is  to  force  the  flow  to  attain  the  correct  streamlines  In  this  vicinity. 

Loads  on  these  elements  are  not.  included  in  the  drag  and  lift  total* 
for  the  entry  body. 

No  load  elements  are  placed  on  a surface  extending  from  the  position  on 
the  model  where  separation  occurs  to  the  effective  planar  surface.  The  actual 
location  of  each  element  is  adjusted  in  successive  runs  until  the  pressures  on 
it  Is  at  a desired  level  (C  ■ 0 for  u vented  cavity).  Typical  results  for 
a disk  cylinder  are  shown  lR  Figure  13  assuming  a vented  cavity.  In  this  case 
the  cavity  was  modeled  using  a single  ring  of  elements  extending  from  the  edge 
of  the  face  to  the  effective  planar  surface.  Although  the  procedure  for  locating 
the  no  loads  elements  is  not  automated  and  therefore  somewhat  tedious,  the 
results  do  account  for  much  of  the  difference  between  experiment  and  theory. 
Fortunately,  it  is  generally  not  necessary  to  include  the  cavity  at  timeB  prior 
to  the  peak  load. 

Correcting  Pressure  on  Modified  KlementH 


During  the  vertical  entry  of  axisymmetric  bodies  it  is  appropriate  to 
apply  a correction  to  the  pruaHurea  calculated  on  modified  elements.  Under 
these  conditions  the  body  is  grldded  with  elements  having  a pair  of  sides 
parallel  to  the  water  surface.  Following  the  step  size  rule  of  equation  (24) 
elements  are  submerged  In  exactly  two  steps.  On  odd  numbered  steps  the  elements 
adjacent  to  the  water  surface  are  nil  modified  while  in  even  numbered  ones  they 
are  not  modified,  Apparently,  the  pressure  levels  predicted  in  the  odd  steps 
are  not  consistent  with  those  calculated  in  the  even  ones.  TIiIb  is  illustrated 
by  considering  the  vertical  entry  of  a 22,5,  45  and  70  degree  half-angle  cones.. 

The  drag  coefficients,  non-dimens ional ized  by  the  local  surface  diameter,  are 
given  in  Table  1 as  a function  of  depth.  Since  this  problem  is  conical  in  nature 
the  drag  should  be  the  same  at  each  depth.  For  the  first  few  steps  error  may  be 
expected  since  the  entire  cone  is  being  modeled  with  only  a few  elements.  However, 
results  should  converge  to  a common  value..  It  is  clear  from  Table  1 tliut  the  odd 
and  even  numher  steps  are  converging  to  different  values,.  In  order  to  determine 
the  better  of  the  two  answers,  calculated  pressure  distributions  are  compared 
to  experiment  for  typical  odd  and  oven  numbered  steps  in  Figures  14  and  15., 

These  pressure  distributions  are  very  similar  except,  near  the  water  surface. 

Figure  14  clearly  shows  that  the  pressure  on  the  modified  element  Is  too  large,. 

A simple  correction  factor,  F ,,  can  be  determined  which  when  multiplied  by  the 
pressure  on  the  modified  element  brings  the  total  drag  calculated  in  even  and 
odd  steps  into  lino  with  one  another.  The  value  of  this  correction  fuctor  has  been 
plotted  in  Figure  16  as  a function  of  depth  for  the  three  different  cones  under 
consideration. 

When  applying  the  present  model  to  the  vertical  entry  of  axisymmetric  bodieB, 
either  the  drug  calculated  in  odd  numbered  steps  should  be  ignored  or  the  correction 
factor  F _ should  be  used.  In  the  remainder  of  this  report  a value  of  .67  is  used. 

The  above  problem  does  not  arise  during  oblique  entry  since  the  edges  of  the 

generated  elements  are  not  parallel  with  the  water  surface  and  the  number  of 

modified  elements  Is  fairly  constant  from  step  to  step.  This  does  not  mean  that 

silt'll  a correction  Is  not  necessary.  However,  there  Is  no  systematic  method  tor 

chosing  F(.  on  general  bodies,  Experience  suggests  that  on  slender  bodies  (nose 

length/D  >l)n  value  of  ,67  should  he  used  and  In  all  other  oblique  entry  cases 

F should  lie  set  to  unitv.. 
c 
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APPLICATION  OF  THE  CODE  TO  SPECIFIC  EXAMPLES 

The  previously  described  computer  program  has  been  applied  to  the  oblique 
entry  of  disk  cylinders,  ogives,  spheres,  and  to  the  vertical  entry  of  ogives, 
cusps,  cones,  and  spheres.  In  this  section  these  calculations  are  compared 
to  experimental  results.  In  assessing  the  validity  of  the  water-entry  model 
it  should  be  kept  in  mind  that  some  uncertainty  exists  with  regard  to  many 
of  the  experimentally  determined  quantities.  Also,  measured  quantities  may 
not  be  equivalent  to  calculated  ones.  The  measured  pressure  represents  the 
value  at  a specific  point  on  the  model  while  the  calculated  result  reflects  the 
average  for  an  element  of  finite  size.  These  two  quantities  become  synonymous 
on  elements  well  below  the  water  surface. 

Vertical  Entry  of  Axlsymmetric  Bodies 

In  this  section  the  vertical  entry  of  cones,  cusps,  spheres,  and  ogives 
is  considered.  The  two-dimensional  nature  of  thesn  problems  insures  that  a 
airuile  value  of  C,  accurately  characterizes  the  rate  of  surface  wetting  about 
the^entire  perife^y  of  the  model.  Also,  separation  of  flow  on  cones  and 
cusps  can  be  categorically  ruled  out  until  after  the  shoulder  of  the  model 
has  entered  the  water.  These  cases  thus  provide  an  ideal  opportunity  for 
testing  the  proposed  predictive  method. 

Vertical  cone  entry  calculations  are  compared  to  the  experimental  results 
of  Baldwin^, 14  which  were  taken  at  entry  velocities  of  1C  to  32  ft/a ec. 

Using  a correlation  developed  in  thlB  work,  an  expression  for  the  wetting 
factor  can  be  obtained  which  is  applicable  to  cones  with  a half  angle  greater 
than  7.5  degrees! 


C 5 T (25) 

w [1.  - ,3960  + . 2079  - .1240  J] 

c c c 

Here  0 is  the  cone  half  angle  in  radians.  Measured  and  calculated  pressure 
coefficients  on  22,5,  45,  and  70  degree  half-angle  cones  are  shown  In  Figures 
15,  17,  and  18.  The  experimental  values  represent  a correlation  based  on 
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Baldwin,  J.  L.,  An  Experimental  Investigation  of  Water  Entry,  PhD  Dissertation, 
U.  of  Maryland,  1972 

^baldwin,  J„  L. , "Vertical  Water  Entry  of  Cones,"  Naval  Surface  Weapons  Center, 
White  Oak  Laboratory,  NOLTR  71-25  (1971) 
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conical  similarity  which  has  been  corrected  to  reflect  a constant  entry  velocity 
Use  of  normalized  depth  as  the  independent  variable  is  appropriate  since  the 
computational  model  also  produces  conically  similar  results.  Excellent  agreement 
is  obtained  between  calculations  and  experiment,.  In  particular,  predicted 
pressure  coefficients  reflect  the  reversal  in  functional  form  exhibited  by  the 
data  with  increasing  cone  angle..  Experimental  data  near  the  tip  of  the  cone  is 
not  Bhown  since  Baldwin  has  indicated  that  there  were  an  absence  of  measurements 
in  the  region^. 

In  Figure  19,  the  calculated  drag  of  finite  length  cones  ure  compared  to 
Baldwin's  results*  flood  agreement  exists  up  to  the  point  where  the  cone  becomes 
completely  submerged  which  coincides  with  the  occurrence  of  peak  drag.  At 
later  times  the  calculated  values  are  too  low.  Improved  agreement  between 
experiment  and  theory  would  probably  be  obtained  if  the  cavity  were  modeled. 

The  calculated  drag  on  vertically  entering  cusps  and  ogives  are  compared 
to  the  experimental  measurements  of  reference  (16)  in  Figures  20  to  23,  The 
dimensions  of  these  bodies  are  given  in  their  respective  figures.  Predicted 
values  are  in  good  agreement  with  experiment  prior  to  the  drag  peak,,  To 
accurately  determine  the  drag  peak  on  the  cusp  models  the  grid  was  extended 
past  the  actual  shoulder  by  one  row  of  elements*  The  present  calculative 
method  anticipates  the  end  of  the  cusp  one  atep  before  it  occurs  making  this 
procedure  necessary  (see  Equ,  (18)),  At  times  following  the  point  of  peak 
drag,  forces  on  the  entry  body  are  calculated  both  with  and  without  a simulated 
cavity*  In  the  cases  depicted  in  Figures  20,  21,  and  23,  inclusion  of  the 
cavity  brings  the  computed  drag  into  close  agreement  with  experiment*  The 
formation  of  a cavity  does  not  appear  important  for  the  ogive  of  Figure  22, 

This  body  is  the  slenderest  of  the  four  models  with  little  surface  discontinuity 
at  the  shoulder-afterbody  Junction, 

A systematic  method  for  determinine  C in  the  above  four  cases  involves 
substituting  the  local  body  angle,  f)  5,,  in¥o  equation  (25)  to  determine  C aB 
u function  of  time..  It  would  seem  plausible  to  use  either  the?  local  ungYe 
on  the  effective  planar  surface  or  at  the  original  surface.  The  validity  of 
these  two  approaches  can  he  examined  for  the  cusp  modelH  whore  the  peak  drag 
can  be  assumed  to  occur  as  the  shoulder  of  the  model  is  wetted  Tile  better 
of  these  two  methods  Is  the  latter,  but  even  it  overpredicts  Cw  resulting  in  a 
premature  wetting  of  the  shoulder,  A correction  factor  can  be  determined  which 
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produces  the  actual  time  of  shoulder  wetting  when  multiplied  by  the  factors 
calculated  from  equation  (25)  using  0£  at  the  original  surface.  The  resultB 
shown  in  Figures  20  and  21  reflect  correction  factor  values  of  .97  and  .94 
respectively.  Correction  factor  with  values  greater  than  unity  might  be 
postulated  for  ogives  since  their  profiles  are  convex  instead  of  concave. 

However,  this  type  of  adjustment  was  not  carried  out  and  the  computations 
for  the  ogive  models  used  the  values  of  Cy  defined  by  equation  (25).  Examination 
of  Figures  22  and  23  indicates  that  such  a correction  would  have  reduced  the 
discrepancy  between  theory  and  experiment. 

Calculations  for  the  vertical  entry  of  a sphere  are  compared  to  Hisewanger's^ 
experimental  measurements.  These  tests  were  made  at  an  entry  velocity  of  23.5  ft/ 
sec.  Pressures  were  measured  at  a number  of  points  on  the  model  surface  and 
integrated  to  produce  total  drag.  The  response  times  of  successive  gages  were 
used  to  give  the  following  expression  for  the  wetting  factor: 

C - 1.  736  - . 829/t*-  (26) 

w 

Calculated  pressures  are  compared  to  measured  ones  in  Figures  24  through  26.  Only 
the  pressure  measurements  made  while  the  transducers  were  fully  wetted  are  shown. 

The  predicted  stagnation  pressure  is  over  estimated  at  early  times  but  in  good 
agreement  otherwise.  At  intermediate  distances  from  the  stagnation  point  the 
predicted  pressure  is  below  the  measured  one.  However,  far  from  the  stagnation 
point,  as  is  shown  in  Figure  26,  the  calculated  pressure  is  again  close  to 
experiment.  The  predicted  drag  is  compared  to  experiment  in  Figure  27  with  beBt 
agreement  being  obtained  at  early  times.  The  measured  drag  does  not  account  for  ^g 
model  deceleration.  However,  these  results  are  in  good  agreement  with  Mosteller's 
constant  velocity  data. 

Oblique  Entry  of  Arbitrary  Bodies 

The  oblique  entry  of  arbitrary  bodies  constitutes  a more  rigorous  test  of 
the  predictive  method  since  these  cases  are  three  dimensional,.  Calculations 
are  compared  to  experiment  for  cones,  disk  cylinders,  spheres,  and  ogives. 


^Niaewanger , C,  R.,  "Experimental  Determination  of  Pressure  Distribution  on  a 
Sphere  During  Water  Entry",  NAVWEPS  Report  7808,  Oct  1961 
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Predicted  and  measured  pressure  distribution  on  a 45  decree  half-angle  cone 
entering  vertically  but  at  an  agle  of  attack  are  shown  in  Figure  28.  The 
experimental  data  are  unpublished  results  of  .1.  L.  Baldwin  of  NSWC/WOL  taken 
on  the  windward  and  leeward  ray  of  the  cone  at  incidences  of  10  and  20  degrees. 
Equation  (25)  was  used  to  determine  C . On  the  windward  ray  the  cone  angle 
was  incremented  by  the  angle  of  attach  while  on  the  leeward  side  it  was  decreased 
by  this  amount.  As  can  be  seen  from  this  figure,  results  are  generally  in  good 
agreement  with  measurements. 

Experimental  data  on  the  oblique  entry  of  disk  cylinders  can  be  found  in 
the  work  of  Norman-^,  Hosteller^®,  and  Baldwin^ representing  entry  velocities 
of  25  to  325  ft/sec.  Based  on  the  latter  two  sources  and  data  taken  in  the 
present  study,  a C value  of  1.45  is  selected  for  all  entry  angles.  Existing 
Information  for  this  parameter,  shown  in  Figure  29,  contains  extensive 
scatter  and  hence  this  choice  is  a rough  estimate.  In  Figure  12,  the  calculated 
drag  is  compared  to  Baldwin's  empirical  correlation  of  experimental  data.  Both 
theory  and  experiment  agree  quite  well  over  a wide  range  in  entry  angle. 

Calculated  pressures  are  compared  to  Aronson's  experimental  data  In  Figure  30 *• 

In  Figures31and  32  calculations  are  compared  to  Aronson's  pressure  measurements 
on  a three-inch-diameter  ogive  cylinder  entering  obliquely  at  100  ft/sec  and  an  angle 
of  60  degrees.  This  body  has  a flat  face,  1,5  inches  in  diameter,  and  rounded 
shoulders  with  radii  of  .75  Indies.  The  mean  measure  value  of  1.36  is  used  in  the 
computations.  Experimental  results  Indicate  that  pressures  rise  more  quickly, 
to  a higher  peak  and  fall  more  rapidly  on  the  lower  portion  of  the  model  face. 

This  is  particularly  evident  for  measurements  made  on  the  shoulder  of  the  model. 

The  computed  results  closely  reflect  this  change  in  the  pressure  trace  associated 
with  transducer  location.  However,  peak  pressures  are  consistently  overpredicted. 
Fortunately,  the  high  peak  pressures  act  on  extremely  small  areas  and  thuB  have 
little  effect  on  the  actual  load  for  the  size  elements  being  used.  To  illustrate 
this  point,  experimental  and  calculated  loads  on  the  computational  element  nearest 
to  the  center  of  the  face  are  plotted  in  Figure  33.  The  experimental  load  Is 
obtained  by  applying  the  data  from  location  1 in  Figure  32.  The  pressure-time 
history  of  each  point  in  the  computational  element  is  assumed  described  by  this 
relation. 
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A brief  summary  of  this  experimental  work  will  soon  be  available  from  the 
National  Technical  Information  Service  in  a report  titled  "Prediction  of 
Surface  Pressures  During  Water  Impact"  by  Wardlaw  and  Aronson 


23 


/ 


1 1-  I l-il.  tVjil. 


it  I II  I •hiiTliifll 


NSWC/WOL/TR  77-16 


In  Figure  27  the  calculated  drag  for  the  oblique  entry  of  a sphere  is  plotted 
against  experimentally  smoothed  curves  given  by  May  which  are  constructed  from 
the  data  of  references  (18),  (20),  (21),  and  (22).  This  information  reflects 
entry  velocities  between  11  and  225  ft/sec.  The  wetting  factor  for  the  oblique 
entry  of  a sphere  has  not  been  extensively  investigated.  A value  of  1.35  was 
selected  based  on  White 's23  limited  results.  Reasonably  good  agreement  is 
obtained  between  calculated  and  experimental  values. 

Calculated  loads  on  a slender  ogive  body  entering  obliquely  are  compared 
in  Figures  34  and  35  to  the  unpublished  drag,  normal  force  and  pitching  moment 
data  by  Baldwin.  The  wetting  factor  for  this  case  does  not  appear  to  have  been 
investigated  experimentally.  A value  of  1.1  is  used  which  corresponded  to  the 
average  value  obtained  for  the  vertical  entry  of  this  body  using  equation  (25). 

Best  agreement  between  measurements  and  calculations  occurred  using  the 
OGIVE  grid  option,  discussed  in  Appendix  C,  which  produces  elements  with  a pair 
of  edges  parallel  to  the  water  surface.  Consistent  with  previous  discussion, 
the  calculated  loads  at  odd  numbered  steps  are  discarded.  Analytical  results  are 
in  closest  agreement  with  experiment  at  the  entry  angle  of  75  degrees  as  shown 
in  Figure  35.  The  premature  decrease  in  the  calculated  drag  at  6 ■ 45°  in 
Figure  34  can  probably  be  attributed  to  the  formation  of  a cavity  along  the 
upper  surface  of  the  ogive.  The  use  of  no  load  elements  to  model  the  water- 
cavity  Interface  could  conceivably  decrease  this  discrepancy.  Consistent 
with  trends  visible  in  previous  examples,  the  underprediction  of  the  drag  initially 
occurs  near  the  point  where  the  total  load  on  the  slender  ogive  model  peaks. 

The  present  calculative  procedure  does  not  include  the  contribution  to  normal 
force  and  pitching  from  the  formation  of  an  underpressure  cavity.  Hence  in  cases 
where  this  effect  is  Important  the  calculated  normal  force  and  pitching  moment 
will  not  be  very  accurate. 
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SUMMARY  AND  CONCLUSIONS 

This  technical  report  outlines  a systematic  method  for  calculating  surface 
pressures,  forces  and  moments  on  arbitrary  bodies  during  the  early  phases  of 
water  entry,  A potential  flow  model  is  assumed  and  the  free  surface  is  approximated 
by  an  effective  planar  surface  empirically  located  at  the  splash  height.  The 
computational  techniques  of  Smith  and  Hess  are  used  to  solve  the  the  potential 
problem.  This  requires  that  the  surface  of  the  body  be  described  by  planar 
quadrilateral  elements.  Using  Bernoulli's  equation  the  average  pressure  is 
calculated  on  each  element  and  then  integrated  to  produce  total  forces  and  moments 
on  the  entry  body.  Through  the  use  of  no  load  elements  it  is  possible  to  model 
the  cavity  which  form  about  the  entry  body,  but  this  is  generally  not  necessary. 

The  described  method  of  calculation  has  been  applied  to  a number  of  different 
cases  in  which  experimental  data  is  available.  For  vertical  entry  this  includes 
cones  with  and  without  angle  of  attack,  ogives,  cusps,  and  spheres.  The  oblique 
entry  case  has  been  studied  for  disk  cylinders,  spheres,  blunt  and  slender 
ogives  covering  entry  angles  between  30  and  90  degrees.  The  predicted  pressure 
traces  accurately  duplicate  experimental  results,  reflecting  not  only  overall 
body  geometry  but  also  location  on  the  body  surface.  The  calculated  loads  are 
in  good  agreement  with  experimental  values,  particularly  prior  to  the  point  of 
peak  loading.  At  later  times  no  load  elements  must  be  used  to  model  the  water- 
cavity  interface. 

Although  the  current  predictive  method  is  a viable  engineering  tool,  some 
shortcomings  are  evident,  Most  notably,  pressures  on  elements  adjacent  to  the 
water  surface  are  often  overpredicted.  This  is  not  surprising  considering 
the  singularity  which  exists  at  the  water  surface  in  the  current-  formulation. 

An  empirical  correlation  scheme  based  on  an  experimental  data  correlation  (e.g,, 
reference  (24))  might  offer  substantial  improvement.  Finally,  it  is  clear  from 
the  studied  examples  that  positive  steps  must  be  taken  to  model  the  cavity  after 
the  point  of  peak  load.  Provisions  should  be  made  for  automating  this  procedure. 


^Baldwin,  J.  L. , and  Steves,  H.  K.,  "Vertical 
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TABLE  1 CALCULATE  CONE  DRAG  AS  A FUNCTION  OF  STEP  NO. 
AT  ODD  NUMBER  STEPS  ELEMENT  ADJACENT  TO  THE 
WATER  SURFACE  IS  MODIFIED  WHILE  ON  EVEN 
NUMBER  STEPS  IT  IS  NOT. 
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FIO.  2 ELEMENTS  WITH  SIMILAR  SOURCE  STRENGTHS.  ONLY  HALF  OF  A BODY  WITH  v’  • *' 
PLANE  18  GRIDDED.  EACH  ELEMENT  HAS  A CORRESPONDING  SYMMETRIC,  IMAGE 
AND  IMAGE  SYMMETRIC  ELEMENT  OF  THE  SAME  SOURCE  STRENGTH  MAGNITUDE. 
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FIG.  3 COMPUTATIONAL  GRID.  THE  MODEL  SURFACE  IS  DIVIDED  INTO  PLANAR 
QUADRILATERAL  ELEMENTS.  ALSO  SHOWN  IS  THE  INTERSECTION  OF  THE 
WATER  SURFACE  AND  THE  MODEL  DURING  THE  FIRST  THREE  STEPS. 


NSWC/WOL/TR  77-16 


3.0  r- 


0.6 


FIG.  7 DRAG  OF  A DISK  AT  AN  ENTRY  ANGLE  OF  60  DEGREES  AS  A FUNCTION  OF  Cw 
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FIG.  13  PROFILE  OF  THE  CAVITY  ABOUT  A OISK  CYLINDER  A1  SEVERAL  ENTRY  ANGLES 

CALCULATED  USING  NO  LOAD  ELEMENTS. EFFECTIVE  PLANAR  SURFACE 

---WATER-CAVITY  INTERFACE. 
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FIG.  16  PRESSURE  DISTRIBUTION  ON  A 46  DEGREF  HALF  ANGLE  CONE  ENTERING  VERTICALLY, 
THE  SHADED  CIRCLES  REPRESENT  THE  CALCULATED  VALUES  AT  ELEMENT  CENTROIDS 
WHILE  THE  HORIZONTAL  LINES  INDICATE  THE  EXTEND  OF  EACH  ELEMENT,  THE  SOLID 
CURVE  IS  EXPERIMENTAL  DA  TA  BY  BALDWIN13  7HE  ELEMENT  ADJACENT  TO  THE 


WATER  SURFACE  IS  NOT  MODIFIED.  Cw  « 1 46. 
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FIG.  17  PRESSURE  DISTRIBUTION  ON  A 22.6  DEGREE  HALF-ANGLE  CONE  ENTERING  VERTICALLY 

THE  SHADED  CIRCLES  REPRESENT  THE  CALCULATED  VALUE  AT  EACH  ELEMENT  CENTROID. 
WHILE  THE  HORIZONTAL  LINE8  INDICATE  THE  EXTENT  OF  EACH  ELEMENT.  THE  SOLID 
CURVE  18  EXPERIMENTAL  DATA  BY  BALDWIN13,  Cw  - 1.14. 


43 


4 


WATER  SURFACE 


NSWC/WOL/TR  77-16 


FIG.  18  PRESSURE  DISTRIBUTION  ON  A 70  DEGREE  HALF-ANGLE  CONE  ENTERING  VERTICALLY. THE 

8HADED  CIRCLES  REPRESENT  THE  CALCULATED  VALUE  AT  EACH  ELEMENT  CENTROID  WHILE  THE 
HORIZONTAL  LINES  INDICATE  THE  EXTENT  OF  EACH  ELEMENT.  THE  SOLID  CURVE  IS 
EXPERIMENTAL  DATA  BY  BALDWIN13,  C - 1.38. 
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L 20  CALCULATED  AND  MEASURED  DRAG  ON  A VERTICALLY  ENTERING 
MEASURED16  ^CALCULATED  WITH  A CAVITY  SIMULATED  BY 
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FIG.  30  CALCULATED  AND  MEASURED  PRESSURE  TIME  HISTORIES  AT  TWO  DIFFERENT 
POSITIONS  ON  THE  SURFACE  OF  A DISK  CYLINDER,  0-80  AND  Vf»  100  FT/SEC. 
— — MEASURED  AT  POSITION  1 <r  - 0,098  p - 4°) •CALCULATED  AT  POSITION 
1.--  -MEASURED  AT  POSH  ION  2 <r  - O'  ^CALCULATED  AT  POSITION  2. 
MEASUREMENTS  BY  ARONSON. 
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FIG.  32  CALCULATED  AND  MEASURED  PRESSURE  TIME  HISTORIES  AT  THREE  DIFFERENT  POSITIONS 
ON  THE  8URFACE  OF  AN  OGIVE  CYLINDER, ©60  AND  Vj  - 100  FT/8EC.  — MEASURED  AT 
POSITION  1 (r  - 0)  A CALCULATED  AT  POSITION  1.— MEASURED  AT  POSITION  2 (r-O.MB',0-  90  ) 
• CALCULATED  AT  POSITION  2.-*-MEA8URED  AT  POSITION  3<r  - 0.112’,  0-  90  ) ■ CALCULATED 
AT  POSITION  3.  DATA  ARE  BY  ARONSON 
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FIG.  33  CALCULATED  AND  EXPERIMENTAL  LOAD  ON  THE  CENTER  ELEMENT 

OF  THE  OGIVE  CYLINDER  MODEL.  A CALCULATED.  —EXPERIMENTAL. 
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APPENDIX  A 


FORMULAS  FOR  THE  INFLUENCE  OF  PLANAR,  QUADRILATERAL  ELEMENTS 

This  section  lists  formulas  for  calculating  the  Influence  of  planar  elements 
on  arbitrary  polntaln  space,  x,y,z,  These  expressions  are  taken  from  reference  (11) 
and  are  listed  here  for  the  sake  of  completeness.  The  symbols  to  be  used  are 
defined  In  Figure  A-l»  Calculations  are  carried  out  using  the  element  coordinate 
system  and  hence  the  velocity  components  are  referenced  to  these  exist 

V5  “ "S12Q12  “ S23Q23  " S34Q34  “ S41Q41 

Vn  " C12Q12  + C23Q23  + C34Q34  + C41Q41 


whore  t 


sign(z)  [A0  — Jj^  ” ^23  "*  ^34  ” ^41^ 
$12  + $23  + $34  + $41  " I * I Ae 


0 unless  R^2*R23  ,R34  *R41  > 0 Then  A4>  Is  2ir 
RijQi.)  + |a|jij 


- - (y  - VC.J 


n + rj  - dld 


si  s1 

agn(Rlj)[tan"1(|i2-  \-&  ).  - tan"1*^  |-ii  )] 
R1J  r1  Rij  ri 

[(x  - q)2  + (y  - n^2  + z2)l/2 


- x)cij  + (nk  - yjs^ 

5,  - 5, 


n,  - n. 


KC,  - q>2  + (n  i - n,)2]1/2 


NSWC/WOL/TR  77-16 


Ac  large  distances  from  the  element  (r,5/t  > A),  the  element  may  be  treated 
as  a point  source,.  This  greatly  simplifies  calculation  of  the  velocity  and 
potential: 


T 

♦ " ro 


v ■ z 

2 V 

1 • T <C3  - «l><r|2  - V 


Vr,  • fj  1 
o 

The  equations  for  $ and  V can  bn  evaluated  in  any  coordinate  system.  The  final 
velocity  components  will  be  referenced  to  whatever  Bystem  is  used.  At  intermediate 
distances  from  the  element,  2,45  < r /t  < 4,  mult ipole  expansions  are  used  to 
evaluate  the  influence  coefficients,0 


+ “ I00w  + 2(I20\x  + 2Illwxy  + I02wyy) 
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APPENDIX  B 

DESCRIPTION  CT  THE  COMPUTER  PROGRAM 

The  computer  program  consists  of  the  main  routine,  ENTRY,  eighteen  subroutines 
and  one  function.  The  flow  chart  for  the  program  is  given  in  Table  B-l  and  the 
main  program  variables  are  listed  in  Table  B-2.  A program  listing  is  provided 
in  Table  B-3. 

Program  execution  is  initiated  in  ENTRY  by  calling  subroutine  INDATA  which 
reads  the  required  input  (see  Appendix  C for  a description  of  input  cards). 

Based  on  options  specified  by  the  user  a grid  is  set  up  on  the  surface  of  the 
model  using  subroutines  STANG,  LISTG,  and  OGTVG.  Node  points  describing  the 
surface  of  the  body  are  stored  in  body  (x*,  y’,  z')  and  water  surface  coordinates 
(x,  y,  z)  in  arrays  X and  XT  respectively.  The  identification  numbers  of  the 
four  nodes  defining  the  Jth  element  are  stored  in  array  IN(1,  J),  IN(2,  J),  IN(3,  J) 
and  IN(4,  J).  The  coordinates  of  each  node  point,  nodes  making  up  each  element 
and  element  centroids  are  printed  in  subroutine  ELEMST,  This  terminates  the 
element  definition  portion  of  the  program  which  is  executed  only  once. 

The  following  sequence  is  executed  at  each  step.  The  center  of  gravity  of 
the  entry  body  is  inserted  the  prescribed  increment  in  depth  and  the  corresponding 
position  of  each  node  point  is  calculated  by  subroutine  ADVN,  Subroutine  CREL 
is  called  which  examines  all  elements  to  determine  which  are  submerged.  Elements 
split  by  the  water  surface  are  redefined  using  new  node  points  which  located  on 
the  water  surface  and  stored  at  the  ends  of  arrays  X and  XT.  The  nodes  making  up 
submerged  elements  are  stored  sequentially  in  array  IT,  The  original  reference 
number  of  each  submerged  element  is  stored  in  array  IQ  and  the  code  indicating 
its  state  (i.e,,  modified,  unmodified  or  split)  in  array  IM.  Upon  completion 
of  subroutine  CREL,  MAN  is  called  which  calculates  the  area,  centroid  and  unit 
vectors  of  each  element.  If  the  PRINT  option  is  used,  element  information  is 
printed  out  by  ELEMST  for  elements  which  are  modified  or  unmodified  for  the 
first  time.  This  information  includes  centroid  and  node  locations  in  both 
coordinate  systems,  element  -ireas  and  unit  vector  components.  The  arrays  X,  XT, 

IN,  XC,  XCP,  VNO,  E,  IT,  PHIO,  IQ,  IM  are  written  on  TAPE  15  by  ENTRY.  The 

coefficients  of  matrices  A,  B,  C,  and  D are  calculated  by  subroutine  AliAT  using 

the  equations  of  Appendix  A and  written  on  TAPE  16,  To  determine  the  induced 
velocity  normal  to  each  element  centroid  (i.e,,  the  right  hand  side  of  equation 
(13)) subroutine  VNORM  is  called.  The  source  strengths  of  equation  (13)  are 

calculated  using  DDECOMP  and  DSOLVE,  These  subroutines  are  called  from  SUAS 

which  sets  up  matrix  A in  appropriate  blocks  using  information  on  Tape  16. 

The  element  information  stored  on  TAPE  15  and  the  matrices  B,  C,  and  D Btored 
on  TAPE  16  are  recalled  to  allow  velocity  and  potential  values  to  be  calculated 
at  each  centroid  element  in  subroutine  OUTPUT  using  equation  (14),  If  the 
PRINT  option  is  used,  this  information  is  also  printed.  Data  necessary  for 
calculating  pressures  and  loads  are  stored  on  TAPE  17, 
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The  above  procedure  is  repeated  until  termination  occurs.  This  can  be 
triggered  by  reaching  a step  number  greater  than  IMAX,  having  the  'number  of 
submerged  elements  exceed  NPIM,  or  having  the  computational  time  approach  the 
Job  time  limit.  This  last  option  is  necessary  since  pressures  and  loads,  of 
principal  concern,  are  calculated  at  the  end  of  the  program,  To  insure  that 
these  quantities  ere  computed  the  code  estimates  the  amount  of  time  required 
to  complete  each  atep  before  starting  it.  If  the  total  estimated  time  estimated 
exceeds  90  percent  of  the  job  time  limit,  the  flow  field  calculation  is  terminated 
and  the  program  branches  to  subroutine  PRESF  which  calculates  pressures  and 
loads  using  the  equations  (18)  through  (21), 
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TABLE  B-2 

H 

'.j 

:| 

MAIN  VARIABLES 

- 

VARIABLE 

TYPE 

DEFINITION 

i 

AN,AX,AY, ANS, 
AXS.AYS 

ARRAY 

TEMPORARY  STORAGE  IN  AMAT  AND  OUTPUT 

■ ' d 

CGL 

SIMPLE 

z'  COORDINATE  OF  CENTER  OF  GRAVITY 

CM 

ARRAY 

WETTING  FACTORS  USE  DURING  CONSTANT 
ORIENTATION  ENTRY 

' 1 
■! 

CWT(I) 

ARRAY 

WETTING  FACTOR  BETWEEN  STEPS  1*1  AND  I 

j 

D 

SIMPLE 

MODEL  DIAMETER.  USED  ONLY  IN  CALCULATING 
DIMENSIONLESS  QUANTITIES 

"3 

% 

Hi 

DT(I) 

ARRAY 

INCREMENT  IN  TIME . BETWEEN  STEPS  1-1  AND  I 

• 

i 

E(I, J,K) 

ARRAY 

ELEMENT  UNIT  VECTORS.  I*  1,2,3  ARE 
COMPONENTS  ALONG  x,y,z  AXIS.  J - 1,2,3 
ARE  UNIT  VECTORS  £,  n.  AND  y RESPECTIVELY. 
K IS  ELEMENT  NUMBER 

1 

! 

HMIN 

SIMPLE 

INITIAL  MODEL  DEPTH 

‘‘j 

r 

HHAX 

SIMPLE 

INCREASE  IN  MODEL  DEPTH  DUE  TO  RISING 
MOTION  OF  SURFACE 

# 

ICON 

SIMPLE 

IF  ICON  - 0,  CONSTANT  ORIENTATION  ENTRY 
IS  ASSUMED.  IF  ICON  » 1 VARIABLE 
ORIENTATION  IS  USED 

I EM 

SIMPLE 

NUMBER  OF  SUBMERGED  ELEMENTS 

IM(I) 

ARRAY 

I IS  ELEMENT  NUMBER,  IM  - 0,  IM  - 1,  AND 
IM  - 2,  INDICATE  NOT  MODIFIED,  MODIFIED 
AND  SPLIT  ELEMENTS  RESPECTIVELY 

IN(I,J) 

ARRAY 

NODES  DEFINING  CORNER  LOCATIONS.  I IS 
THE  CORNER  NUMBER  AND  J IS  THE  ELEMENT 
NUMBER 

INP 

SIMPLE 

TOTAL  NUMBER  OF  DEFINED  NUDES.  INCLUDES 
NODES  GENERATED  TO  DESCRIBE  MODIFIED 
ELEMENTS 

IP 

ARRAY 

USED  IN'  DIJLCOMP  AND  DSOLVE 

B-4 
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TABLE  3-2  (Continued) 


VARIABLE 

TYPE 

DEFINITION 

IPRINT 

SIMPLE 

IF  IPRINT  - 1 PRINT  OPTION  IP  EXERCISED 

IPHI 

SIMPLE 

IF  IPHI  - 1 PLANAR  SYMMETRY  IS  ASSUMED 

iq(K) 

ARRAY 

ORIGINAL  IDENTIFICATION  NUMBER  OF 
SUBMERGED  ELEMENTS.  K TS  THE  IDENTIFICATION 
NUMBER  OF  THE  SUBMERGED  ELEMENT 

ITCI.K) 

ARRAY 

IDENTIl'I CATION  NUMBER  OF  NODES  DESCRIBING 
THE  Kth  SUBMERGED  ELEMENT 

NCW 

SIMPLE 

NUMBER  OF  WETTING  FACTORS  TO  BE  USED 
DURING  CONSTANT  ORIENTATION  ENTRY 

NDIM 

SIMPLE 

MAXIMUM  ALLOWABLE  NUMBER  OF  SUBMERGED 
ELEMENTS 

NEL 

SIMPLE 

NUMBER  OF  ELEMENTS  INITIALLY  DEFINED 

PHI 

ARRAY 

TEMPORARY  STORAGE  IN  AMAT 

PHIS 

ARRAY 

TEMPORARY  STORAGE  IW  OUTPUT 

SIG(K) 

ARRAY 

■SOURCE  STRENGTH  OF  ELEMENT  K 

STOR 

ARRAY 

TEMPORARY  STORAGE  T.N  SUAS 

SUMT 

SIMPLE 

SO  PERCENT  OF  TIME  LIMIT 

VENTRY 

SIMPLE 

ENTRY  VELOCITY 

VEX(I) ,VEY (I) , VEZ(l) 

APRAY 

VELOCITY  COMPONENTS  IN  THE  x,y,  AND 
z DIRECTION  USED  BETWEEN  STEPS  I nnci  1-1 

VNO(I) 

ARRAY 

VELOCITY  NORMAL  TO  ELEMENT  I 

WX(I) 

ARRAY 

ANGULAR  VELOCITY  IN  THE  PITCH  PLANE 
BETWEEN  STEPS  T AND  1-1 

X(I,K) 

ARRAY 

COORDINATES  OF  THE  Kth  NODE.  I - 1,2,3 
REFER  TO  THE  x\y\z'  AXES  RESPECTIVELY 

XC( I ,K) 

ARRAY 

COORDINATES  OP  THE  CENTROID  OF  THE  Kth 

ELEMENT.  I «■  1,2,3  REFER  TO  THE  x*,y',  AND 
z * AXES  RESPECTIVELY 
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VARIABLE 

XCP(I,K) 


},  J 

K 


I) 


XCPB(I) 

XT(I,IC) 
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TABLE  B-2  (Continued) 


TYPE  DEFINITION 

ARRAY  COORDINATES  OF  THE  CENTROID  OF  THE  Kth 

ELEMENT.  I - 1,2,3  REFER  TO  THE  x,y,  AND 
z AXIS  RESPECTIVELY 

ARRAY  COORDINATES  OF  THE  CENTER  OF  GRAVITY. 

I - 1,2,3  REFER  TO  THE  x,y,z  AXES 
RESPECTIVELY 

ARRAY  COORDINATES  OF  THE  Kth  NODE.  I - 1,2,3 

REFER  TO  THE  x,y,z  AXES 
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APPENDIX  C 

USER  INSTRUCTIONS  AND  SAMPLE  RUNS 

The  current  version  of  the  code  can  be  applied  to  arbitrary  bodies.  It  was 
developed  on  a CDC  6500  and  requires  about  105K  storage  octal.  The  grid 
describing  the  entry  body  may  contain  up  to  750  nodes*  and  500  elements,  however, 
execution  will  terminate  when  more  than  300  of  these  elements  become  submerged. 

This  appendix  described  the  available  program  options,  necessary  input  cards 

and  output  format.  Sample  runs  are  provided  to  illustrate  the  use  of  this  program. 

Program  Options  and  Required  Input 

Program  input  can  be  divided  Into  three  parts.  In  the  first,  the  basic 
program  options  are  specified! 

Card  No.  Variable  Format 

1 CONSTANT  or  VARIABLE  body  orientation  2A4 

2 PRINT  or  DON'T  PRINT  3A4 

3 ASYMMETRIC  or  SYMMETRIC  mode  3A4 

Under  the  CONSTANT  body  orientation  option  the  entry  model  is  assumed  to  retain 
its  initial  orientation  and  velocity  throughout  the  entry  process.  As  discussed 
previously,  the  natural  problem  variable  in  this  case  is  depth  rather  than  time. 
With  little  increase  in  computational  time,  pressures  and  forces  can  be  evaluated 
for  a number  of  different  wetting  factors,  C . The  VARIABLE  body  orientation 
option  allows  the  velocity,  orientation,  wetting  factor  and  time  increments 
between  steps  to  be  varied.  The  only  restriction  iB  that  the  angular  velocity  of 
the  body  must  be  small  enough  to  insure  that  the  depth  of  the  body  increases 
monotonically  in  time.  The  maximum  number  of  steps  is  limited  to  49. 

The  PRINT  option  is  used  to  obtain  flow  field  and  element  information  at  each 
step  of  the  calculation.  It  is  applied  only  for  diagnostic  purposes.  The 
second  option,  DON'T  PRINT,  is  recommended  and  produces  only  grid  information 
and  the  final  pressures  and  forces  on  the  model. 

If  the  SYMMETRIC  mode  option  is  used,  the  entry  model  is  assumed  to  possess 
planar  symmetry  about  the  y-z  plane.  The  ASYMMETRIC  option  does  not  assume  any 
symmetry  and  hence  can  be  applied  to  arbitrary  bodies.  This  mode  is  also  used 
on  symmetric  bodies  where  Vx  is  none-zero. 


*Storage  is  set  up  for  800  nodes.  However,  this  must  also  include  room  for 
nodes  generated  at  each  step. 
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The  second  set  of  input  cards  describes  the  entry  conditions  and  the  required 
information  differs  depending  on  whether  the  CONSTANT  or  VARIABLE  body 
orientation  option  la  used.  For  the  CONSTANT  option  the  following  data  cards 
are  required: 


Card  No. 

4 

5 

6 

7 


Variables 

IMAX,  D,  VENTRY,  ANG,  SUMT,  HMIN,  DH,  ALPHA 

CGL,  FCF,  ANGB,  NNLD 

NCW,  CW(l),  CW(2) (CW(NCW) 

omit 


Format 

15,  5X,  7F10.0 
3F10.0,  15 
15,  5X,  7F10.0/ 
(8F10.0) 


These  variables  are  defined  as  follows: 

IMAX  Number  of  steps  at  which  pressures  and  loads  are  calculated. 

The  present  calculative  procedure  Inserts  the  model  into  the 
water  in  a series  of  steps,  each  at  a greater  depth  than  the 
preceding  one.  When  the  step  count  becomes  greater  than  .IMAX 
execution  is  terminated 

D Diameter  (in  feet).  This  quantity  is  only  used  for  calculating 

force  coefficients 

VENTRY  Entry  velocity  in  ft/ sec 

ANG  Orientation  of  the  model  (in  degrees)  relative  to  the  water 

surface  (see  Fig.  C-l) 

SUMT  Program  time  limit.  This  must  correspond  to  the  time  limit  on 

the  job  card 

HMIN  Initial  body  depth  (i.e.,  measured  from  the  lowest  point  on  the 

body).  This  parameter  is  zero  if  the  loads  are  calculated  from 
time  of  initial  wetting.  Note  that  if  this  variable  is  not 
zero  pressures  and  forces  are  first  calculated  at  HMIN  + 2DH. 

This  parameter  allows  pressures  and  forces  at  a particular  depth 
tc  be  determined  without  calculating  the  entire  force-time  history 
from  Initial  wetting. 

DH  Increment  in  depth  in  feet  between  successive  steps.  It  is 

necessary  to  coordinate  this  variable  with  the  specified  model  grid 
which  is  defined  on  the  last  set  of  data  cards.  The  following 
apply  to  determining  DH: 

a.  OBLIQUE  ENTRY  WITH  STANDARD  GRID  OPTION.  DH  should  be  picked 
so  that  the  average  element  is  submerged  in  two  steps.  On 
models  of  complex  shape  this  criteria  can  only  be  satisfied  in 
the  mean  and  primary  consideration  should  be  given  to  the  portion 
of  the  body  which  experiences  the  greatest  load.  Generally  this 
will  be  on  elements  whose  plane  is  perpendicular  to  the  direction 
of  motion. 
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ALPHA 

CGL 

FCF 

ANGB 

NNLD 

CW 


b.  VERTICAL  ENTRY  WITH  STANDARD  GRID  OPTION  OR  OBLIQUE  ENTRY 
WITH  THE  OGIVE  OPTION,  For  vertical  entry  or  if  the  OGIVE 
grid  option  is  used,  elements  will  have  a pair  of  side  parallel 
to  the  water  surface.  In  this  case  it  is  important  to  choose 
the  step  size  very  precisely  so  that  each  element  will  be 
submerged  in  exactly  two  steps.  To  insure  that  the  top  row 
of  elements  is  included  in  an  unmodified  state  in  the  code  and  that 
the  next  row  of  elements  is  excluded,  the  actual  water  surface 
should  fall  a small  distance  e above  the  upper  edge  of  the  top  row 
of  element  to  be  Included  as  shown  in  Figure  C-2,  Here  0 < e < Ah 
where  Ah  is  defined  by 

Ah  ■ /average'  element  area 
1000 

Angle  of  attack  in  degrees  (see  Fig,  C-l) 

z'  coordinate  of  the  center  of  gravity  (see  Fig,  C-l) 

Pressure  correction  factor  on  elements  with  a modification 
code  of  1.  For  the  oblique  entry  of  blunt  bodies  (nose  length/ 
diameter  < 1)  set  to  unit.  For  other  cases  use  a value  of  .67, 

Yaw  angle  in  degrees  of  Vj.  Velocity  components  in  the  x,y,z  direction 
are  VTsin  (ANGB),  -V.cos(ANG)cos(ANG  + ALPHA)  and  - VTcos(ANG) 
ain(ARG  + ALPHA) 

Number  of  wetting  factors  to  be  used.  Since  the  most  appropriate 
value  may  not  be  clear,  for  little  extra  computational  cost 
pressure  and  loads  may  be  calculated  for  several  different 
wetting  factor  values 

Wetting  Factor.  This  parameter  describes  the  rate  of  surface 
rise  and  is  equal  to  the  ratio  of  h/h'  defined  in  Figure  1. 

For  best  results,  the  test  cases  reported  on  should  be  used 
as  a guide.  An  approximate  rule  for  determining  this  parameter 
is  as  follows: 

(1)  POINTED  BODIES  (ALSO  INCLUDES  SLIGHTLY  BLUNTED  ONES).  Determine 
the  angle,  9 , between  the  tangent  to  the  body  surface  and 
the  body  axis  at  both  the  nosetip  and  base  of  the  nose.  At 
the  nose  neglect  any  effect  due  to  body  blunting.  Insert 
the  two  resulting  values  of  6C  in  radians  into: 

C 5 r (C-l) 

w [1  - .3969  + ,2879  1 - .1240 

c c c 

Average  the  two  calculated  values  of  C to  obtain  the  final 

one  to  be  used  in  the  code.  If  ALPHA  is  non-zero  increment 

9 by  ALPHA 
c 
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(2)  FLAT  PLATES.  Use  a value  of  1.45  for  ANG  > 45  degrees 
and  1.55  for  ANG  < 45. 


(3)  SPHERICAL  BODIES.  Use  a value  of  1,55  for  near  vertical 
entry  and  1.35  for  oblique  entry. 


The  classification  of  an  arbitrary  body  into  one  of  the  above 
categories  is  a matter  of  experience.  On  complex  shapes 
classification  should  be  based  on  the  portion  of  the  body 
sustaining  the  majority  of  the  impact  loading. 


If  the  VARIABLE  body  orientation  option  is  used  the  following  data  cards 
are  required: 

Card  No.  Variables  Format 


4 IMAX,  D,  VENTRY,  ANG,  SUMT,  HMIN 

5 CGL,  FCF,  NNLD 

6 NVP 

7.1  VX(1),  VY(1) , VZ(1),  CW(1),  DH(1) 


15,  5X,  7F10.0 
2F10.0,  10X,  15 

6F10.0 


7.NVr  VX(NVP),  VY (NVP) , VZ<NVP),  WX(NVP),  CW(NVP) , 6F10.0 

DH(NVP) 

The  variables  on  cards  4 and  5 are  defined  above,  In  this  case,  VENTRY  is  only 
used  in  determining  the  force  and  pressure  coefficients  and  ANG  is  the  initial 
body  orientation,  The  body  velocity,  wetting  factor,  and  increment  in  depth  for 
each  step  is  defined  in  cards  7. 


NVP  Number  of  different  Bteps  at  which  entry  conditions  are 

specified 


VX(I),  VY(I),  Velocity  components  in  the  x,y,z  directions  of  the  center  of 
VZ(I)  gravity  in  ft/sec  applied  between  steps  1-1,  and  I 

WX(I)  Angular  velocity  in  degrees/sec  in  the  pitch  (y-z)  plane 

applied  between  stcpB  1-1  and  I 

CW(I)  Wetting  factor  applied  between  the  1-1  and  I step.  If  the 

value  of  this  parameter  remains  constant  from  step  to  step 
use  the  Instruction  for  determining  this  variable  given  in 
the  CONSTANT  orientation  section.  For  the  vertical  entry 
(VX'VY-0)  of  pointed  bodies  an  estimate  of  this  parameter 
for  each  step  can  be  obtained  by: 

a.  Determining  the  depth  of  the  entry  body,  H,  below  the 
original  surface  at  the  Btart  of  the  step 

b.  Calculating  the  angle,  O’-  between  a tangent  to  the  body 
surface  and  the  body  axis,  z'  ■ H 
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c.  Substituting  Into  equation  (C-l)  to  determine  C . Where 

w 

e - 0 ’ + 90  - ANG 
c c 

d,  For  blunt  bodies,  [(nose  length) /diameter]  < .75,  increase 
this  angle  by  7%  on  ogives  end  decrease  It  by  the  same  amount 
on  cusps 

DH(I)  Increase  in  depth  in  feet  of  the  center  of  gravity  between  steps 

1-1  and  I.  See  instruction  in  the  CONSTANT  orientation  entry 
flection 

The  entry  velocity  at  step  I is  taken  to  be  the  average  of  that  at  steps  1-1 
and  I.  It  is  only  necessary  to  specify  data  cards  for  the  first  few  steps  in 
which  the  above  parameters  change.  For  steps  larger  than  NVP  the  parameter 
values  at  step  NVF  are  used. 

The  final  set  of  data  cards  1b  used  to  define  the  grid  on  the  surface  of 
the  entry  body.  The  three  available  options  for  constructing  a grid  on  the 
body  surface  are  STANDARD,  OGIVE  and  LIST.  These  can  be  used  singularly, 
in  combination  with  one  another  and  can  be  called  in  arbitrary  sequence.  The 
only  restriction  la  that  the  lowest  point  on  the  body  should  occur  on  that 
part  of  the  grid  constructed  by  the  first  option  called.  To  indicate  the 
desired  options,  the  following  input  cards  are  required: 


Card  No. 

Variable 

Format 

8 

N 

15 

9.1 

u 

option  1 

3A4 

9!n 

option  N 

3A4 

Here  N is  the  number  of  options  to  be  used.  Recommended  option  are  provided  in 
Table  C-l. 

The  grid  representing  the  surface  of  the  entry  body  should  cover  only  the  nose 
of  the  model  and  not  the  afterbody.  In  all  cases  the  pressures  on  the  afterbody 
are  small.  Furthermore,  on  bodies  with  sharp  shoulders  such  as  a disk  cylinder, 
the  flow  separates  at  the  edge  of  the  model  face.  If  the  afterbody  is  gridded, 
the  flow  is  required  not  to  separate  since  the  invicid  boundary  conditions 
are  enforced  at  the  centroid  of  each  element.  This  is  physically  unrealistic 
and  hence  neglecting  the  afterbody  is  appropriate, 

A description  of  the  three  available  options  follows.  Under  no  circumstances 
should  right  angles  be  modeled  directly.  If  the  body  under  consideration  has 
such  a surface  discontinuity,  it  should  be  modeled  with  a 89.9  or  90,1  degree 
angle. 
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STANDARD 

i This  opclon  is  applicable  Co  axisynunecric  bodies  or  axisymraetric  portions  of 

>J  arbitrary  bodies.  The  user  specifies  rings  along  which  nodes  are  located. 

\ Adjacent  nodes  are  combined  to  form  elements.  A typical  grid  for  a flat, 

i circular  plate  is  shown  in  the  Figure  C-3.  The  required  input  is: 

j Card  No.  Variables  Format 

10  NROWS,  IANG,  I SUP 

11.1  R(l),  2(1) 

11. NROWS  R (NROWS) , 2(NR0WS),  IW(NROWS) 


IANG  If  IANG  ■ 0,  only  half  of  the  face  is  grldded  as  shown. 

If  IANG  * 1,  the  complete  face  is  gridded 

NROW  Number  of  grid  rings 

ISUP  If  I SUP  ■ 1,  the  stagnation  element  (number  1)  is  removed. 

This  option  is  used  for  running  pointed  objects.  For  such 
bodies,  R(l)  should  be  very  small  (i.e.,  D/1000)  but  must 
be  finite. 

If  ISUP  • 0 this  element  is  included. 

R(I)  Radius  of  ring  I in  body  fixed  coordinates  (x',  y',z')  in  feet. 

2(1)  z*  coordinate  of  ring  is  in  feet. 

IW(I)  number  of  elements  in  the  area  between  ringB  I and  1-1, 

Delete  this  variable  on  card  1.  If  IW  - 0,  elements  are 
automatically  selected  so  that  they  are  approximately  9quare, 


315 

2F10,0,15 

2F10.0,15 


•J 


i 

i 
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LIST 


This  option  requires  that  the  user  input  the  list  of  nodes  and  elements 
to  be  used  in  the  run  and  hence  is  applicable  to  arbitrary  bodies.  The  nodes 
can  be  read  in  any  order,  however,  they  are  numbered  sequentially  for  internal 
use  in  the  code.  Each  element  is  constructed  using  nodes  from  the  input  list. 
The  identification  numbers  of  nodes  defining  the  four  corners  of  each  element 
must  be  read  in  a clockwise  order  with  respect  to  an  observer  on  the  outer 
surface  of  the  element.  The  required  input  cards  are: 


Card  No. 

Variable 

Format 

10 

NP,  NE 

215 

11.1 

• 

x’(D,  y'U),  x'd) 

3F10.0 

ii!np 

x'(NP),  y'(NP),  z* (NP) 

3F10.0 

12.1 

IN(1,1) , IN(2,1),  IN(3,1) . IN (4 , 1) 

415 

12.  NE 

IN(1,NE),  IN(2,NE),  IN(3,NE),  IN(4,NE) 

415 

NP 

number 

of  node  points  to  be  read  in. 

NE 

number 

of  elements  to  be  read. 

x' (I).y'(I)  location  of  the  Ith  rode  in  body  fixed 
z'(I)  coordinates  (x'.y'.z')  in  feet 

TN(l.I), 

IN(2,I), 

IN(3,I) , 

IN (4, I) 


Identification  number  of  the  nodes  defining  the 
four  corners  of  the  element 
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OGIVE 

Thla  subroutine  is  applicable  to  pointed  axisymmetrlc  bodies  entering 
obliquely.  It  grids  the  body  surface  with  elements  having  a pair  of  edges 
parallel  to  the  water  surface  as  shown  in  Figure  C-4  , This  feature  is  particularly 
desirable  for  conical  bodies  where  ANG  1b  less  than  90  degrees.  Under  these 
conditions  the  described  code  produces  conical  results  and  the  complete  pressure- 
time history  can  be  obtained  using  calculated  values  at  a single  depth.  This 
subroutine  allows  elements  to  be  packed  near  the  water  surface  and  minimizes 
the  need  to  use  the  local  similarity  assumption.  Use  of  this  subroutine  is 
also  recommended  for  examples  in  which  normal  force  is  of  importance.  Alignment 
of  elements  with  the  water  surface  will  produce  more  accurate  results  than  the 
STANDARD  grid  option.  If  this  subroutine  is  used  to  grid  non-conical  bodies, 
the  constructed  element  will  be  slightly  non-planar,  however,  this  does  not 
seriously  effect  the  calculations. 


Card  No, 

Variable 

Format 

10 

NROWS , NBODY,  IANG 

315 

11,1 

• 

R(l),  Z(l) 

2F10.0 

11. NBODY 

R (NBODY),  Z (NBODY) 

2F10.0 

12.1 

H (1 ) , B(l) , IW(1) 

2F10.0.15 

12.NR0WS 

H ( NROWS ) , B (NROWS) , IW(NROWS) 

2F10.0.15 

Here  cards  11  are  used  to  specify  the  body  shape  while  cards  12  define  the  grid. 
Up  to  SO  points  may  be  used  to  define  the  body.  Between  these  points  the  body 
is  defined  by  linear  interpolation. 

The  above  variables  ore  defined  ast 

NKOWS  Number  of  element  rows 


NBODY 
I ANG 

H (I ) , Z(I) 


Number  of  points  defining  the  body  geometry 

If  IANG  - 0 only  the  right  half  of  the  body  is  gridded.  If 
lANG  " 1 the  entire  body  is  gridded 

2 2 1/2 

Polar  coordinates  of  body  profile  in  feet,  R ■ (x1  + y*  ) , 

Z - z'.  R(l)  and  Z(l)  must  be  zero 


H(I)  z*  coordinate  of  the  upper  edge  of  the  Ith  element  row  in  the 

y'-z'  plane  In  feet 

B(I)  Orientation  in  degrees  of  the  entry  body  axis  with  reBpect  to  the 

water  surface  at  the  end  of  step  21-  If  the  entry  is  under 
constant  condition  B(I)  is  equal  to  ANG 


IW(I) 


Number  of  elements  in  the  Ith  row.  Must  be  specified. 
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SAMPLE  RUNS 

Four  sample  runs  have  been  provided  to  Illustrate  the  uBe  of  various  code 
options.  In  each  case  a brief  discussion  is  given  of  the  entry  problem  and 
its  particular  peculiarities.  This  is  followed  by  a listing  of  the  input  cards 
and  the  resulting  output  . In  the  final  section  of  this  appendix  a discussion 
is  given  of  the  output  format.  Table  C-2  gives  a brief  outline  of  the  four  cases 
to  be  present. 

Example  1 

A disk  cylinder  entering  at  60  degrees  and  100  ft/sec  is  modeled  using 
a coarse,  12-element  STANDARD  grid.  Consistent  with  the  preceding  discussion, 
the  grid,  shown  in  Sketch  3,  covers  only  the  face  of  the  model.  Noting  that  the 
distance  between  successive  grid  ring  is  .05  ft,  an  increment  in  depth  of 
•0125  feet  is  chosen  insuring  that  the  average  element  1b  submerged  in  two 
steps.  Consistent  with  instructions,  Cw  and  FCF  are  Bet  to  1,45  and  1 respectively. 

Example  2 

The  vertical  entry  of  a 45-degree  half-angle  cone  at  20  degrees  incidence 
is  studied  using  an  OGIVE  grid  with  10  rows  of  elements  and  8 elements  In  each 
row.  This  problem  is  conically  similar  indicating  that  pressures  and  loads  are 
adequately  defined  by  considering  only  a single  depth,  hence  IMAX  ■ 1.  HMIN  is 
selected  as  .47814  ft  which  places  the  water  level  Just  above  the  8th  element 
row  and  DH  is  set  at  .02988  ft  to  insure  that  each  element  is  submerged  in  two 
steps  exactly.  Pressures  are  calculated  at  HMIN  + 2DH  which  corresponds  to  a 
level  just  above  the  9th  element  row.  Two  values  of  C are  used.  These  are 
arrived  at  by  substituting  0 ± ALPHA  in  equation  (C-l);  The  larger  value  is 

most  appropriate  on  the  windward  ray  of  the  cone  while  the  smaller  applies  to 
the  leeward  ray. 

Example  3 

The  vertical  entry  of  an  ogive  traveling  at  50  ft/sec  is  considered  under 
the  variable  entry  option.  A STANDARD  grid  consisting  of  90  elements  is  con- 
structed on  the  body  nose.  The  entry  velocity  and  orientation  are  fixed  through- 
out but  the  wetting  factor  is  allowed  to  vary  from  step  to  step.  The  value  of 
this  parameter  is  determined  using  the  procedure  outlined  for  calculating  CW 
under  the  VARIABLE  entry  option.  The  increments  in  depth  between  successive 
steps,  DH ( I ) , are  selected  with  care  to  insure  that  the  water  surface  is  a small 
distance,  e,  above  the  upper  edge  of  elements  adjacent  to  the  water  surface  at 
even  numbered  stepB.  In  examining  the  total  calculated  drag  a higher  degi.ee 
of  reliability  is  placed  on  result  obtained  at  even  numbered  steps. 
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PROGRAM  OUTPUT 

Under  the  Initial  heading  of  I,X, Y,Z,XC,YC,ZC  are  listed  all  nodes., 

Column  1 la  the  node  Identification  number,  columns  X,Y,Z,  are  the  node  positions 
in  the  x'.y’,*'  coordinates  and  XP,  YP,  ZP  are  the  node  pooitlons  in  x,y,z 
coordinates  when  h • 0,  The  second  section  of  the  output  lists  the  constructed 
elements.  The  integers  under  the  heading  ELEMENT  are  the  element  reference 
numbers,  XC,YC,ZC  are  the  coordinates  of  the  element  centroid  referenced  to  the 
axes  x*,y’,z’  and  RC  and  TC  are  the  radial  location  measured  from  the  z*  axis 
and  the  angular  position  measured  from  the  +y'axis  of  the  element  centroid,.  The 
third  section  of  output  indicates  numbers  of  the  completed  steps  and  the  time 
required  to  execute  that  particular  step.  If  a problem  occurs  during  execution 
and  the  program  does  not  terminate  normally,  the  step  in  which  the  problem 
occurred  can  thus  be  determined. 

The  remainder  of  the  output  defines  the  loading  on  body  at  each  step  in  the  entry 
of  the  model.  At  each  step  the  program  calculates  the  time  and  dimensionless 
time  (Vjt/D)  from  initial  impact.  This  is  followed  by  a listing  of  the 
submerged  elements.  The  reference  number  (i.e.,  ref,  no,)  of  an  element  is  the 
original  identification  number  assigned  to  it  while  its  number  (no.)  is  a 
temporary  identification  parameter  for  the  particular  depth  in  question.  The 
integer  under  the  heading  of  MOD  identifies  whether  the  element  iB  modified  or 
split.  The  codes  0,1  and  2 refers  to  unmodified,  modified,  and  split  elements 
respectively.  The  area  of  each  element  Is  that  of  its  submerged  portion.  Also 
listed  for  each  element  are  its  pressure  coefficient (CP) , centroid  depth  in 
dimensionless  time  (T*) , pressure  (P),  and  total  load  (force).  The  output  for  each 


depth  is 

concluded  by  providings 

FX 

force  along  the  x axis  2 

CX 

X force  coefficient  (FX/((tt^“)  1/2  p V^) 

FD 

drag  force 

CD 

drag  force  coefficient  **FD/((ttD^/4)-  p V^) 

FN 

normal  force  (i„e„,  in  the  -y*  direction) 

FD 

2 

normal  force  coefficient  (FN/((~-)l/2  p 

) 

SMX 

x component  of  moment  taken  about  the  model 
x moment  coefficient  SMX/ ((^jp  )|  p V^) 

center  of 

gravity 

MX 

SMY 

y component  of  the  moment  about  model  center 

of  gravity 

MY 

■nn3  i 2 

y morent  coefficient  SMY/((— jr— )tj-  P ) 

SMZ 

z component  of  moment  taken  about  the  model 

center  of 

gravity 

MZ 

irn3  i ? 

z moment  coefficient  UMZ/((-^-)y  p ) 
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r 


TYPE  OF  BODY 

Axlaymmetrlc 

Pointed 

Axlaymmetrlc 

Blunt 

All  other  bodies 


TABLE  C-l 

RECOMMENDED  GRID  OPTIONS 

RECOMMENDED  OPTION  COMMENTS 

Vertical  Entry  Oblique  Entry 

STANDARD  OGIVE 

OGIVE 

STANDARD  STANDARD 


LIST  LIST  On  pointed  bodies, 

it  is  advisable  to 
construct  elements 
with  edges  parallel 
to  water  surface  as 
done  in  OGIVE 


SAMPLE  RUNS 
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EFFECTIVE 

PLANNAR 

SURFACE 


FIG.  C-2  PROFILE  OF  CONE  GRID 
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FIG.  C-3  GRID  OF  A CIRCULAR  PLATE 


OAT*  CAROS  FOB  EXAMPLE 
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OAT*  CAROS  FOR  £*A*Pt_E 
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